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Abstract. We study the non-equilibrium dynamics of the spherical fcrromagnet 
quenched to its critical temperature, as a function of the magnetization of the initial 
state. The two limits of unmagnetized and fully magnetized initial conditions can be 
understood as corresponding to times that are respectively much shorter and much 
longer than a magnetization timescale, as in a recent field theoretical analysis of the n- 
vector model. We calculate exactly the crossover functions interpolating between these 
two limits, for the magnetization correlator and response and the resulting fluctuation- 
dissipation ratio (FDR). For d > A our results match those obtained recently from 
a Gaussian field theory. For d < 4, non-Gaussian fuctuations arising from the 
spherical constraint need to be accounted for. We extend our framework from the fully 
magnetized case to achieve this, providing an exact solution for the relevant integral 
kernel. The resulting crossover behaviour is very rich, with the asymptotic FDR X°° 
depending non-monotonically on the scaled age of the system. This is traced back 
to non-monotonicitics of the two-time correlator, themselves the consequence of large 
magnetization fiuctuations on the crossover timescale. We correct a trivial error in our 
earlier calculation for fully magnetized initial states; the corrected FDR is consistent 
with renormalization group expansions to first order in 4 — d for the longitudinal 
fluctuations of the 0{n) model in the limit n — > oo. 



1. Introduction 

The use of fluctuation-dissipation ratios (FDR) has proved very fruitful in the last 
decade or so for quantifying the non-equilibrium dynamics of glasses and other systems 
exhibiting aging. In the context of mean-field spin glass models with infinite-range 
interactions, the FDR, commonly denoted X, has been used to formulate a generalized 
fluctuation-dissipation theorem (FDT) where X is interpreted in terms of an effective 
temperature, Tg^ = T/X for the slow, non-equilibrated modes of the system [I]. The 
properties of X and Tefj have attracted much attention, based on the hope they might 
allow a generalized statistical mechanical description for a broad class of non-equilibrium 
phenomena [21 [3l H] . 
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However, the generalized FDT can be shown to hold exactly only for infinite- 
range models. A matter of recent intense interest has been whether the appealing 
features of this mean-field scenario survive in more realistic systems with finite-range 
interactions |2]. A class of systems that has proved useful in this context is represented 
by ferromagnets quenched from high temperature to the critical temperature Tc or below 
(see e.g. Refs. [3 [6l O [8], |9] and the recent review [ID]). The non-equilibrium dynamics 
in these systems is due to coarsening, i.e. the growth of domains with the equilibrium 
magnetization (for T < Tc) or equilibrium correlation structure (for T = Tc), and slows 
down as domain sizes increase. In an infinite system, equilibrium is never reached, 
leading to aging; the age- dependence of two-time quantities has a simple physical 
interpretation in terms of the growth of the domain lengthscale [11] . Coarsening systems 
therefore provide a physically intuitive setting for the study of aging phenomena as 
observed e.g. in glasses, polymers and colloids. They are, of course, not completely 
generic; compared to e.g. glasses they lack features such as thermal activation over 
energetic or entropic barriers. 

We focus in this paper on critical coarsening, i.e. coarsening at Tc, where interesting 
connections to dynamical universality exist. The FDR X is determined from correlation 
and response functions which, in aging systems, depend on two times: the age t„ of the 
system and a later measurement time t. In contrast to mean-field spin glasses, where X 
is constant within each "time sector" {e.g. t — t^ = 0{1) vs t — t^ growing with t^), in 
critical coarsening the FDR is a smooth function of t/t^. This makes the interpretation 
of T/X as an effective temperature less obvious. To eliminate the time- dependence one 
can consider the limit of times that are both large and well-separated. This defines an 
asymptotic FDR 



An important property of this quantity is that it should be universal [5l [IDj in the sense 
that its value is the same for different systems falling into the same universality class 
of critical non-equilibrium dynamics. This makes a study of X°° interesting in its own 
right, even without an interpretation in terms of effective temperatures. 

An intriguing theoretical question which has been addressed recently is whether 
different initial conditions can lead to different universality classes of critical coarsening. 
Due to the universality of X°°, these can be uncovered by studying the effect that 
different initial conditions have on the FDR. Of particular interest has been the effect 
of an initial magnetization on the ensuing coarsening. For the Ising model in high 
dimension or with long-range interactions [12J, one finds that magnetized initial states 
do produce a different value of X°°. This suggests a different dynamical universality class 
from conventional coarsening from unmagnetized states, even though the magnetization 
decays to zero at long times. Further steps in this direction were taken in our recent 
calculation of exact FDRs for magnetized coarsening below the upper critical dimension 
in the spherical model [13]. The propagation of a trivial error meant that the results 
were at variance with the renormalization group (RG) result of Ref. [H] derived for the 
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longitudinal fluctuations in the n oo limit of the 0{n) model within an expansion 
around (1 = 4. We give the corrected results in this paper, and these are consistent with 
the RG calculations (see Appendix A ). This suggests that the equivalence between the 
dynamics of the spherical model and the large-n limit of the 0{n) model extends beyond 
the regime of Gaussian fluctuations, where it is trivial to establish. 

Recently it was emphasized in the context of a fleld-theoretic analysis |15] that one 
should think of the nonzero initial magnetization as introducing a new timescale in the 
system. The two limits of unmagnetized and magnetized initial conditions can then 
be understood as corresponding to times that are respectively much shorter and much 
longer than this magnetization timescale, and one can in fact interpolate between these 
two limits using a crossover function that depends on times scaled by the magnetization 
timescale. This crossover function was calculated in [15] in the classical (Gaussian) 
regime, i.e. above the upper critical dimension, but so far there are no predictions 
for this function for lower dimensions where the critical behaviour is governed by non- 
mean-fleld exponents. We provide the flrst results of this kind in this work by calculating 
the relevant crossover functions exactly for the spherical model in 2 < d < 4, for the 
correlator, response and FDR of the magnetization. 

In Sec. [2] we recall the known crossover behaviour of the magnetization (which is 
directly related to a function g(t)) and the general relations encoding the consequences of 
this for the magnetization correlation and response functions. As in [13] non-Gaussian 
spin fluctuations are important and will be accounted for via the kernel L. Key to 
our analysis for the more complicated functions g{t) in our current scenario is an 
exact solution of the integral equation deflning L that applies independently of the 
time regime. In Sec. [3] we then evaluate the magnetization correlator and response 
for d > 4. As expected, we flnd here full agreement with the Gaussian fleld-theoretic 
calculations [15]. Sec. H] deals with the more interesting case d < 4. Here the analysis 
is more complicated but we can still derive exact results for the asymptotic FDR 
The relevant crossover functions display unexpected non-monotonicities that, close to 
the lower critical dimension d = 2, turn into singularities at intermediate values of the 
scaled system age. We study carefully the relevant scaling regimes for ci — > 2, and 
investigate how they arise from the behaviour of the two-time magnetization correlator. 
Our results are summarized in Sec. El 



2. Setup of calculation and exact solution for L^^^ 

We start by recapitulating briefly the relevant elements of our previous analysis of critical 
coarsening in the spherical ferromagnet [l3]. The model consists of N spins Si on a d- 
dimensional cubic lattice, with sites rj and Hamiltonian H = ^ — Sj)"^ [IS]. The 

spins are real-valued but subject to the spherical constraint J2i Sf = N. Langevin 
dynamics leads to a simple equation of motion for the Fourier components = 
J2iSiexp{-iq-ri) of the spins, dtS^ = -(cjq + z(t))5'q + ^q where tUq = 2^f^;L(l-cosg„) 
is abbreviated to u below and ^q is independent Gaussian noise on each wavevector 



Fluctuation- dissipation relations in critical coarsening 



4 



q = {qi,. . . ,qd), with (^q(t)^*(t')) = 2NT5{t - t'). The Lagrange multipher z{t) 
enforces the spherical constraint; as explained in [13j, it is in reality not just a simple 
function of time but a dynamical variable with fluctuations of 0{N~^/'^) that cause all 
the non-trivial effects in the behaviour of global observables. In terms of the function 
g{t) = exp (2 j/;dt'z{t')^ the Fourier mode response is 



di^w) ^-^(t-t^) ^ m{t) ^-^{t~t^) 
\ 9{t) m{t„) ^ ^ 



In the second equality we have used that the time-dependent magnetization can be 
written as m{t) = (l/A^)(5o(t)) = Ro{t,0){l/N){So{0)) = tuq /^fg(t) with mo = 
(l/A^)(S'o(0)) the initial magnetization. The full, unsubtracted two-time correlator 
Cq{t,tvj) = (l/A^)(S'q(t)5'*(tw)) can be related to its equal-time value by the response 
function, 

Cq(^7^w) = R(ii'ti'tw)Cci{tw,t„) (3) 

The relevant equal-time value is given by 

CM, t) = + 2T /* dt' ^e-2-(*-*') (4) 

git) Jo g{t) 

The function g{t) is determined from the spherical constraint, which imposes 
f{dq)Cq{t,t) = 1. Here and below we abbreviate (dq) = (iq/(27r)'^, where the integral 
runs over the first Brillouin zone of the hypercubic lattice, i.e. q G [— tt, tt]'^. The 
resulting integral equation for g{t) is 

git) = |(rfg)Cq(0,0)e-2-* + 2T f^dt'git')fit-t') (5) 

with fit) = Jidq) e^^'^*. Our first task will be to understand how the solution of this 
crosses over between the magnetized and unmagnetized cases. In terms of the Laplace 
transform gis) = dt git)e~'^^, equation ([5]) reads 

9is)-^^fidq)^^ (6) 

We take as the initial condition the standard choice pTl [T5] of a small magnetization 
mo but otherwise uncorrelated spin fluctuations. The initial equal-time Fourier mode 
correlator can then be written as 

Cq(0,0) = 5q,o^m2 + (l-m2) (7) 

This unsubtracted correlator is C(l) for q 7^ but O(A^) for q = 0. (For the fluctuation- 
dissipation behaviour we will need to look at the connected correlator Cq, which is 
discussed below.) Equation ([7]) yields, bearing in mind that the integral idq) is really a 
sum over the N discrete wavevectors with weight each, 

idq)^^ = ^ + il~ml)fis) (8) 
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Using this in one has at criticahty, where T = = [J (dq) 1/uj] 



[2/(0)]- 



with 



K. 



eq 



Tc (dq) 



+ {l-ml)f{s) 



1 



(9) 



(10) 



uj{2uj + s) 

the Laplace transform of the equihbrium form fl38l) of the kernel K defined below. As 
before [13] we want to look at the long-time limit of g{t), corresponding to small s in ([9]). 
In this regime Keq{s) is given ioi d > 4 by -ft'eq(O) — Keq{s) = as^'^~'^^^'^ and for ci < 4 by 
i^eq(s) = bs^'^-^^/\ with a and b some (i-dependent constants. In the remaining square 
bracket of only the first term is present for a fully magnetized initial state (mo = 1); 
conversely, only the second survives for the unmagnetized case (mo = 0). To see the 
crossover between these limits the two terms need to be of the same order. Because we 
are interested in small s and /(O) is nonzero, this implies that and s must be of the 
same order. We then find to leading order in these small quantities 



9{s) 



^eq'(O) 

„(4-d)/2 



s^ 



mn 



+ 



s 

/(O) 



(rf>4) 
(d<4) 



(11) 



or in the time domain 



'cq 



KZ}{Q) mlt + f{Q) 



m. 



0_^{d-2)/2 _^ 



/(O) 



9{t) 

'V{d/2Y ' r((d-2)/2) 

One can combine these two expressions as 



^(d-4)/2 



(rf>4) 
(rf<4) 



where we have defined 



g{t) = —t-\mlt + c) = —t-'' — + 1 



K 




(rf<4) 
(rf>4) 



Aid 



/Ltd 



t 



hT{d/2) 
i^eq(O) 



(rf<4) 
{d>A) 



(12) 



(13) 



(14) 



and c = (1 — fi;)/(0). In the second equality of f|T3l) we have taken out the factor of c to 
identify the crossover timescale 

(15) 



ml 



which as anticipated in the introduction depends on the initial magnetization of the 
system. In the time domain, our statement of the relevant long-time scaling ni^ ~ s 
can now be phrased as follows: we will be considering the limit of large t, t^ and Tm 
(corresponding to small mo) at fixed time ratios m„ = t^/r^^ and Ut = t/r^^^. For ease of 
comparison with the work of [15] we will write simply u = and mostly work with u 
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and the time ratio x = t/t^ = Ut/u instead of u and Ut- In terms of these variables one 
can write the function g{t) as 

g{t) = —r'^i^ux + l) (16) 
For the magnetization one then finds 

Mt) = ^ = ,(^^^ = C>/^ (17) 



g{t) V '^m C ^/uX^^ Y + 1 



with the exponent a defined as a = 1 — k as in [13]. The last square root equals unity 
for long times if the initial magnetization is kept finite and nonzero (so that u ^ 1). 
Otherwise it gives the well-known correction to the fully magnetized result when the 
initial magnetization is small, i.e. when t ~ Tm [17j. In particular, for Ut = ux ^ 1, 
the magnetization displays critical initial slip, increasing as m(t) ~ t'^/^, before crossing 
over to the decay around Ut = l. Our analysis for the fully magnetized case in 
is now recognized as relating to the limit t,t^ ^ t^, and accordingly all results in this 
paper should reduce to the ones in [13] in the limit u —>■ oo. (Loosely speaking, one can 
think of this limit as corresponding to — > 0, i.e. "mo = oo" [15j.) In the opposite 
limit t,tw ^ Tm we should get back the results for the unmagnetized case mo = 0. In 
terms of our scaling variables, this limit corresponds to u at fixed x. Note that there 
is in principle a third, "mixed" regime where the earlier time t„ ^ but the later time 
t ^ Tm, i.e. u 1 and ux 3> 1. We will see, however, that essentially no new behaviour 
arises here and the crossover between the magnetized and unmagnetized cases, which 
the analysis below will allow us to elucidate explicitly, is governed principally by u. 

We next explore how the crossover effects in g{t) modify the expressions for the 
long-time behaviour of the connected Fourier mode correlator C'q(t, t„) = Cq(t, t„) — 
(l/A^)(S'q(t))(S'*(tw)) = Cqit,t„) — N6cifim(t)m(t^) and the response function Rq(t,t^). 
(Here, as previously, we will not write explicitly the dependence on t^.) From [13] 
we know that the equal-time connected correlator has the same expression as the 
unsubtracted correlator 

C'q(tw, tw) = ^ [(5q(0, 0)e-^-*- + 2T, dt' e-2-(*"-*')(7(t')] (18) 

except for the appropriately modified initial condition Cq(0, 0) = 1 — which - in 
contrast to the unsubtracted Cq - is C(l) for all q. For the zero Fourier mode one sees 
that in the long-time limit the first term is subleading and the integral diverges at the 
upper end so that one can use the asymptotics of g(t'), giving 

u/{2-k) + 1/ (1 

fi'(tw) -'0 U + 1 

Similarly in the ratio of nonzero and zero mode correlators, expressed in terms of the 
scaling variable w = ut^, 

C'q(tw,tw) _ {l-ml)e-^^ + 2T,t^J^dze-^^^'-y^gizt, 



C'o(tw,tw) = -^[C7o(0,0) + 2T, dt'git')] = 2T,t^ "/^^ ^) + 1/(1 ^) (^g) 
tw Jo u + 1 



Co(iw, iw) 1 + 2Tctw lo dz g{zt^ 



(20) 
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one can neglect the non-integral terms for long times and gets 
C^{t^,t^) ^ J^dze-''^^'-yh-%zu + l) ^ J^dze-^^('~^h-^{uz + l) 
C'o(tw,tw) J^dzz-{uz + l) n/(2-K) + l/(l-K) 

Putting the last two results together yields the general scaling 

T 2?/; /■! 

C'q(t„,t„) = -J'c{w,u), J^c{w,u) = —— / dze-'^^'-'h-^iuz + l) (22) 
00 u + 1 Jo 

One checks easily that J-'c{w,u) reduces to the analogous scaling functions for the 

unmagnetized and fully magnetized cases [13] in the appropriate limits u and 

u — > oo. The magnetization response function is the zero mode response Rq. From ([2]), 

using the scaling of the magnetization found in (fTTI) . it is given by 



m{t^) \l ux + 1 

The results above are valid within the Gaussian approximation for the spin 
dynamics in the spherical model, where the small fluctuations in the Lagrange multiplier 
z{t) are neglected. As we saw in [13], in order to study the FD behaviour of the 
magnetization [i.e. of the zero Fourier mode, which is a global observable) when 
an initial nonzero magnetization is present, we need to account for non- Gaussian 
corrections arising from these Lagrange multiplier fluctuations. Fortunately our earlier 
expressions [13] for the resulting magnetization correlator and response are valid for 
arbitrary initial conditions and can be used directly. The magnetization correlator 
including non-Gaussian effects is [I3] 

C(t,M = C«(t,tw) + C(2)(t,t„) (24) 

with 

C^'\t,t^) = C'o(t,tw) - / dt' [M{t,t')Coit„,t') + M{t^,t')Coit,t')]m{t') 



and 



+ J dt'dtlM{t,t')M{t^X)Ht')Ht'JCo{t'X) (25) 
= J dt'dt'^[6{t-t') - M{t,t')m{t')] 

X [5(tw-C)-M(t.,C)m(C)]C'o(t',C) (26) 

C^'\t,t^) = \j dt' dtlM{t,t')M{t^,t'JCC{t',t'J (27) 

where CC(t',t'^) = J{dq) C^(t',t'^). The corresponding expression for the global 
magnetization response including non-Gaussian effects is [13] 

R{t, t^) = J dt' [6{t - t') - M{t, t')m{t')]Ro{t', t„) (28) 

The key function M appearing here is defined as follows. One starts from the kernel 

K{t, t„) = I (dq) i?q(t, tw)Cq(t, tw) = / (dq) R^t, t^)C^{t^, t„) (29) 
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and its inverse L defined by 

J dt' K{t, t')L{t', tw) = 5{t - tw) (30) 

The behaviour of K{t, t„) near tw = ^ can be shown to imply the following structure for 
L 

L{t, t„) =S'it- t„) + 2TJit - t„) - (t, t^) (31) 

where the first term arises from the fact that K{t,t^) is causal [i.e. it vanishes for 
tw > t) and has a unit jump at t^ = t. Finally, M is defined to be proportional to the 
integral of L: 

M{t, t„) = m(t) y* dt' L{t', t„) (32) 

In our previous analysis [13] we had found long-time scaling forms of L^"^^ separately 
for the unmagnetized and magnetized cases, with different methods needed for d > 4 
and d < 4. With the function g{t) no longer being a simple power law, it seems 
difficult if not impossible to adapt these methods to our current crossover calculation. 
Fortunately, however, there is a general and fully exact solution for L*^^^ which applies 
in any dimension and for any g{t). To obtain this, we essentially integrate by parts 
in fl30l) . In the derivative of K with respect to the earlier time argument we separate 
off the contribution from the unit step and write 

dt^Kit, t„) = -5{t - t^) + K'{t, t^) (33) 

where K' , like K, vanishes for t^ > t and is finite elsewhere. Correspondingly we split 
off the first term from ( I3TI) and write 

dt' Lit', t^) = 6{t - tw) + N{t, t„) (34) 



where explicitly 

N{t, t„) = 2T, - f dt' (t', t^) (35) 



and N{t, t„) also vanishes for t^ > t. Integrating by parts in (1301) and substituting these 
definitions then yields 

K'{t, t^) + f dt' K'{t, t')N{t', t^) - N{t, t^) = (36) 

J t-w 

The point of this transformation is that non-equilibrium effects manifest themselves 
in K' in a very simple form. To see this, note from (jl]) for the unsubtracted 
correlator that <9t„Cq(tw, ^w) = -[g'it^)/ gi^^) + 2cj]Cq(tw, ^w) + 2Tc, while from 
dt^Rl{t,t^) = [g'{t„)/g{t^)+2uj]Rl{t,t^). Inserting into m gives 

K'it, t^) = 2% I (dq) Rlit, t^) = ^2tJ [dq] e-2-(*-*w) = -^-^K'^^{t - t„) (37) 
J 'I g{t) J g{t) 

where 

i^eq(t-tw)= /(dg) ^e-2-(*-*-) (38) 
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(with Laplace transform given by (ITU]) ) is the equihbrium form of K{t,t^). With the 
simple multiplicative structure of fl37j) one can now solve the integral equation fl36l) for 
by inspection: 

iV(t,M=iVeq(t-t„)^ (39) 

9[t) 

where Neq{t — t^) is the solution of the equilibrium version of ( |36l) . which is related to 
the corresponding by 

The last approximation gives the scalings for large time differences t — t^, derived 
from the corresponding asymptotic behaviour of L^. The latter is L''^(t — t„) = 
Xd{t - t^Y'^^^^/^ ind < 4 and Lg)(t - t^) = Xd{t - t^f'^-'^^/^ in ci > 4, with a d- 
dependent coefficient [I3]. This behaviour can be derived from the Laplace transform 
of L^q) , which from the equilibrium versions of (I30f3ip follows as 

Lf^{s) = s + %-l/k,q{s) (41) 

Note that A^^eq decays to zero ioi d<A because Lg)(0) = /o°°rfrLg)(r) = 2T^ exactly, 
while ioi d> A it approaches the nonzero limit 2Tc — Lf^{Qi) = l/fXd [13j . 
The kernel M is directly related to from (1321) and (1351) : 

M{t, t„) = m{t) [5{t - t„) + N{t, t^)] (42) 

and in our current context we do not then need to compute L^^^ explicitly. Briefly, 
though, the general solution for L^"^^ is 

9{t) ^ 9\t) 
and we outline in Appendix B| how this retrieves all of our previous results in the 



L(2) (t, t„) = - d,N{t, t^) = ^Lg) (t - t^) + ^^'iZ Ne^it - ^w) (43) 



appropriate limits. The key advantage of the above solution method is that it 
automatically accounts for all non-equilibrium effects by reducing the problem to an 
equilibrium calculation at criticality, where all functions depend only on time differences 
and the relevant integral equation can easily be solved by Laplace transform as shown 
in (jlTj) above. 

With the general solution for M{t,t^), and hence for the magnetization correlator 
and response, now in hand we analyse separately the cases d > 4 and d < 4. 

3. Crossover behaviour in > 4 

We first consider the situation d > A above the upper critical dimension. We expect 
to find here the same results for universal quantities as in the Gaussian field theory 
of [15]. The zero Fourier mode Gaussian correlator and response are obtained from ( fT9l) 
and ( l23ll by setting k = 0: 



Ro{t, tw) = J , C'o(tw, tw) = Tc (44) 
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With these, one has from (15^ . (HQ]) and 



M{t,t')m{t') 



6{t - 1') + 



gin 
git) 



1 

jjd d 



2\n 



-{t-t') 



/x(4-d)/2 



t„ (mX + 1)1/2(^1/ + 1)1/2 



t^'K^ - y) 



+ 



uy + 1 

UX + 1 



1 

/Ltd d 



2\d 



4 



(45) 



where we have rescaled the times with t^ and introduced the scahng variable y = t'/t^. 
In the long-time limit the first and the third terms in the above expression are subleading 
for (i > 4 so 

u{uy + iy/^ 



By inserting this expression into ( l28l) one finds the magnetization response 
R{t, tw) = Ro{t, tw) - J dt' M{t, t')m{t')Ro{t', t^) 



(46) 



u+ 1 



1/2 



U[U 



+ 1)1/2 



X - 1] 



r u + 1 

\ux + 1 



3/2 



MX + 1 y {uX + 1)3/2 

The magnetization correlator is found from ( !25l) and reads after rescaling all times 

m(m+ 1)1/2 



(47) 



C«(t,t. 



T f 



u + 2 
u + 1 



dy 



u 



\j UX + 1 

uyiuy + 2; 



dy- 



[UX 



+ 1)3/2 



dy 



uyiuy + 2) 



(m+1)3/2(mx + 1)1/2 



+ dy dy. 



MX +1)3/2(m + 1)1/2 

u^y^iuy^ + 2) 



(MX +1)3/2(^+1)3/2 



+ dy dy 



+ dy dy. 



Tt 



2 + 3m + 2m2 + iM3 



u^y (m?/„ + 2) 

(mX + 1)3/2(^ + 1)3/2 

^^1/(^1/ + 2) 
(MX + 1)3/2(^+1)3/2^ 

1, 



(m + 1)3/2(^0; + 1)3/2 (48) 

The term C*^^) scales as ~ t^~'^^^'^ for 4 < (i < 6, where the integral fl271) that defines 
it can be shown to be dominated by aging timescales, and as ~ t^^ for d > 6, where 
CC{t',t'^) in flTTI) behaves as a short range kernel, so it is always subleading. Thus 
represents the full long-time magnetization correlator for d > 4. 

The t-dependence in the correlator C = C^^^ is the same as in the response 
R and only occurs via the overall factor (mx + l)~3/2 = (^^^ _|_ i)-3/2_ therefore 
cancels in the resulting FDR which follows after a few lines (using dt^ [tv/F{x, u)] = 
{1 + u du — X dx)F{x, u) to calculate dt^C) as 

T,R{t,t^) 4 (m+1)^ 



X(t,t^ 



X^iu) 



(49) 



5t„C(t,t„) 5(m + 1)4 + -, 
Thus, for d > 4 the FDR is t-independent and hence identical to the asymptotic 
FDR X°°(m) = \imt^t^=uTm^i X{t,tw)- It interpolates between 1/2 (for m <^ 1) 
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Figure 1. Normalized magnetization FD plot for dimensionality d above 4, showing 
the normalized susceptibility x versus the normalized correlation C, for different fixed 
values of Ut = ux as indicated in the figure. For ut = Q the plot is a straight line with 
(negative) slope 1/2, as expected from the unmagnetized limit. As Ut is increased the 
initial slope of the plot converges quickly to 4/5, corresponding to the fully magnetized 
limit, and the crossover to the unmagnetized regime occurs at larger time differences 
and eventually becomes invisible on the scale of the plot. 



and 4/5 (for m ^ 1), reproducing in these limits our previous results for the FDRs 
for unmagnetized and fully magnetized initial conditions [I3]. As expected from the 
universality of X°°, our result for the entire crossover function also exactly agrees 
with that calculated from a Gaussian field theory [15]. The FD plot is obtained by 
graphing the normalized susceptibility ^w) = Tcxit, tw) /C(t, t) versus the normalized 
correlator C(t,t„) = C(t,t^)/C(t,t) at fixed Ut = ux and using x (or u) as the curve 
parameter. The factor of Tc is included in the definition of x to make the equilibrium 
FD plot a line of (negative) slope 1. The susceptibility is obtained from R by integration 
as usual, 



(with some obvious abuse of notation in the representation as an integral over u'). The 
results, displayed in Fig. [H show that for ut = 0, the curve is a straight line with 
(negative) slope 1/2, as expected from the unmagnetized limit. As Ut is increased, 
the initial slope of the plot converges quickly to 4/5, which is expected in the fully 
magnetized regime, and the unmagnetized regime gets progressively squeezed into the 
top left corner of the plot where it eventually becomes invisible. Intuitively, this 
is because for ut = t/r^^ ^ 1 we need to move to relatively much earlier times 
u = ^w/tiii ~ 1 in order for the dynamics to be sensitive to the fact that the initial 
magnetization was small. 



rt + rut 

X{t,t^)= dt'R{t,t') = — duR{ut,u') 



Jt„ U J u 



2t^ (ux + iy/^-ju + iy/^ 

5 u {ux + 1)3/2 



(50) 
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4. Crossover behaviour in d < 4 

In d < A the analysis is somewhat more awkward and leads to highly non-trivial 
magnetization FD behaviour as we will see. As before we will find that all relevant 
quantities vary on aging timescales ~ t^ and so we will exploit the relevant asymptotics 
for large time differences throughout. 

One starts by working out the combination M{t,t^)m{t^) appearing in the 
definition of C and R: 

M{t, tw)m(t^) = N{t, t^)m{t^)m{t) = - t^Y'-'^' (51) 

^9{t)g{t^) 9{t) 4 - ci 

_ 2X^t_^^u + l)^ _ _ 1 

- (4-rf)ct-3./2(^3, + i)3/2i^ 1) - ^^m{x,u) (52j 

where we have defined 

Tm{x, u) = ^-ux^^-''^'^ ^ 53 

2 [ux + ly/^ \ X / 

and used 2Ad/id = {4: — d){d — 2)/2 [13j. Equation (153|1 represents the generalization to 
finite u of the scaling function J^m (x) determined in [13] and reduces to the latter in the 
limit M — ^ cxo as it should. Note that in (ISTll we have directly neglected the contribution 
S{t - t^)m{t^)m{t) ~ S{x - l)t-^t-"/2t-"/2 ^ = f-d/2 because it is subleading for 
long times compared to the main 1/t ~ 1 /tw term in (152|) . 
We note briefiy the explicit expression 

that recurs in a number of calculations below; JFa^(1/?/„, u?/w) is obtained from this by 
replacing x ^ 1, y ^ y^. It will also be useful for later to have the asymptotics of 
Tm{x-,u) for large x, which will give the asymptotic behaviour of the correlator and 
thus of 

•Fm(x, u) = ^x^'-'^' ( ^) (55) 



2 \ u 

For the response function the Gaussian contribution is from (!23l) . after setting 
K = (4 - rf)/2, 



Mt^t^) = x^'-''^l\ (56) 
\ ux + 1 



The overall magnetization response is then found simply by inserting (15311 into ( ]28il and 
rescaling the times as before 



I u+l 
ux + 1 



ux /X - 1\ 

ux + 1 \ X 



(58) 
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For u ^ 1, one retrieves the fully magnetized limit calculated previously [13]. On 
the other hand, the unmagnetized limit, m -C 1, coincides with the Gaussian response. 
This is consistent with the fact that non-Gaussian effects in the FD behaviour of the 
(global) magnetization only need to be accounted for in the case of an initial nonzero 
magnetization [13]. The large-x behaviour of i?(t,tw) for general u, which will provide 
the asymptotic FDR, is easily extracted from fl58|) as 



R(t, t^ 



X 



-(d+2)/4, 



'd 



1 

U 



u 



(59) 



We next turn to the correlator. Rescaling the times with t^ in equation fl26|) one 
has for the first part 

C«(t,t^) = r dy f dy^ \5{x - y) - l^A/(f,My) 



5(1 - y. 



l/w 



(60) 



wy w ) 

From ([3]), f[T9l) and fl56|) the Gaussian factor can for long times be written as 

Co{t^y, t^y„) = 2Tct^[f{y, y^)6{y - y„) + f{y^, y)0{y^ - y)] (61) 

= 2%t^[f{y,y^) + e{y^ - !/)(/(?/«, y) - /(y, ?/«))] (62) 

where we have explicitly accounted for the ordering of the time arguments; the 
dependence on y and y^ is through the function 



y_ 

,yw 



k/2 



uy^/{2 - k) + 1/(1 - k) uy„ + l 



uy^ + 1 
{2uy^)/d + 2/{d-2) 

^ (u2/„ + l)V2(«2/ + 1)1/2 

Inserting ( l62l) into (1601) . one can rewrite C^^-* as 



uy 



(63) 

(64) 



C«(M. 

2T t 



dy / dy^ 



6{x-y) - -Tui^^uy) 

Ob 



5(1 - y. 



M 



.uy^ 



X fiy,yw. 
"1 

dy^ 



dy-J^uijy^uy^ [^(1 - y„) - uy^ 

X [/(yw,y) - /(2/,2/w)] (65) 

The decomposition fl62|) is the analogue of the cancellation trick used for the fully 
magnetized case There the analogue of the first line in flB^ vanished identically. 
This is not the case here, but the procedure remains useful because it makes it easier 
to extract the large-a; limit: the first integral (denoted F below) factorizes, and in the 
second one (denoted S) both integration variables y, y^ are < 1 and so <^ a; for large x. 
Using the factorization, the first double integral can be worked out explicitly for generic 
x: 



^= J^^y S{x -y) -^J^Ai{^,uy) 



dyw [5(1 - l/w) - J^M^-^.uy 



y 



(4-d)/4 



X 



{uy + 1)1/2 

,„ {2uy^)/d + 2/{d-2) 
(ny„ + 1)1/2 



(66) 
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^(4-d)/4 

{ux + ly/ 

rl 



1 - 



d-2 



-X 



(4-d)/2_ 



u 



2 ux + 1 Jo 

X / dy„ 5(1 - y^) - ^Af(^,M2/w) 

-'0 

2u/d + 2/{d-2) 



dy{x-yY'-'y' 



,j, {2uy^)/d + 2/{d-2) 
+ 1)1/2 



(67) 



^{4-d)/4 
{ux + 1)3/2 



(« + 1)1/2 



U 



2 (U+ 1)3/2 

The remaining integral produces Beta functions so that 



X 



(4-d)/4 (2u/d + 2/{d-2) 



u T\d/2) 



u + 



2{d-l) 
d-2 



(68) 



(69) 



[ux + 1)3/2 I + 1)1/2 + 1)3/2 r(^) 

The large-x behaviour is obtained by replacing the prefactor with x~^'^~^'^^^^ /u^^'^. The 
second double integral in ( l65l) can be written explicitly as 



S 



1 



X Jo 

1 



dyTM[^,uy) 

guy ^ 2 
\ d 



d-2 



(m?/+ l)l/2(y+ 1)1/2 
_,,(4-rf)/4 , _ 



d d-2 



+ - / dy^J^M[^,uy- 
X Jo 

X dyJ^M[^,uy 

f dy^TAii^ 
X Jo 



y 



(4-d)/4 



y 



{uy^ + 1)1/2 

(i/4 



2My ^ 2 



(u?/ + 1)1/2 ^ ^ 'd-2 

d/4 

, uy^ ' 



S/w 



dyTAiij.uy 



y 



{uy^+iy^ 

(4-d)/4 



2m?/„ _^ 2 



d-2 



d-2 



^' ^/ (m|/ + 1)1/2 

^,3-3^/4^ ^1 



(70) 



2 (m+ 1)1/2(mX + 1)3/2 7o 



rf^ (a;_^)(^-4)/2^('^-4)/4 



V d 



d-2 



y 



2u 2 

+ 



d d-2 



+ 



'd-2' 



^3-3d/4^2 



(m + 1)3/2(uX + 1)3/2 



dy. (1 - y^f'-'^'^ / dy (x - ^)('^-^)/2^('^-2)/2 
dyAi-y.r-'"'y':~'"'i^^+ ' ' 



^2m?/ ^ 2 



V d ' d-2 

dy{x-yr~'y' 



(71) 



V d ' d-2. 

If we can now take the large-x limit, where (x — ?/)('^~'i)/2 ^ x'^"'^^^/^, the integrals can 
be carried out and the prefactor simplifies, giving after a little algebra: 



S 



2,-(2+d)/4 
m3/2(« + 1)3/2 



d-2 

d 



m(m+ 1) 



d 



d + 2 



M+ 1 



d-2r2(d/2) 
^ f(d)" 



.M^ + 1 



(72) 



Gathering the contributions from (169|1 and (1721) , one gets for the large-x limit of the first 
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contribution to the correlator 



2T t 



with 



A 



B 



C 



D 



d-2 
d+2 

2d 
d + 2 

d 

d-2 
2 

d-2 



^3/2(^+1)3/2 

d-2 r\d/2) 
~ r{d) 
d r\d/2) 
2 r{d) 

2{d-l) T\d/2) 



[Au^ + Bu'^ + Cu + D 



d 



(73) 

(74) 
(75) 
(76) 
(77) 



The second contribution to the magnetization correlation comes from ( 127|) . and to 
make progress here we need the long time behaviour of the two-time function CC for the 
current case d < A. Proceeding as for the fully magnetized scenario [13], we find first the 



scaling of the equal-time value CC{t^, t^) = J (dq) C^(t„, t„) = J (dq) u!~'^J^c{^t„, u). 
This is dominated by small cu, where {dq) = adduj uj'^'''~'^^/'^ with the surface area of a 
unit sphere in d dimensions. Rescaling to w = ut^ gives 

CC{t^,t^) = 7d(«)^^^)/^ 7d(«) = J dww^'-'^l'Tl{w,u) (78) 

Normalizing CC{t',t'^) with the equal time value CC{t',t'), one obtains for t' > t'^ in 
terms of the scaling variables y = t' /t^ and = t'^/t^ (by rescaling in the numerator 
to w = ojt'^ and in the denominator to w = ut') 

CC{t',t'J _ g{t'Jt'j,''-''y^ J dww^'^-^y^e-^''^y/y-~''^J^l,{w,uy^) 



CC{f,t') g{t') f(^-'^y^ Jdww(-^-^y^J^l,{w,uy) 
_uy„ + ljdw ?/;('^-6)/2e-2«'fe/2^w-i)jr2 ^y^- 
uy + 1 J dw ty('^~6)/2jF^(w, uy) 

For t' < t'^, on the other hand, one has 
^^(^',0 CC{t'^,t') y^^-'^l^ ^,{uy^) 



CC{t',t') 
So overall 



where 



ccif^^t'jy^^-^y' id{uy) 

CC{t'X 



CC{t',t') 



Q 



(79) 
(80) 

(81) 
(82) 



uy^ 



(j[w,uy„)e 



-2(y/yw-l)«i' 



uy + 1 

ld{uy^) 



J dw w^'^~^'>/'^J^Q{w,uy) 

(d-4)/2 

Q{yw/y,uy) 



for y/y^>l 
for y^/y <l 



(83) 



Idiuy) \y^J 

In the limit u ^ 1, this function should match with Q{x) defined in [13] for fully 
magnetized intial conditions, for x = y/y^. Unfortunately there was a typographical 
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error in the definition of Q{x) as given in [13], wliicli propagated tlirough tlie remainder 
of the calculation. In Appendix A we state the correct versions of all the relevant 
equations. These include, in particular, the first-order expansions around d = 4 and 
d = 2 of X°° in the fully magnetized case. 

Having clarified the scaling behaviour of CC{t',t'^) in d < A, we can now work out 
(7(2) f];om fl271) by multiplying and dividing by m{t')m{t'^) and using (|52|) 



C(2)(t,t^) = 1 f dt' r dt'^——^, :FM{-,uy)TM{- 

^ ' 2 Jo Jo ^ m{t')m{t'Jtt^ 



xCC{t',t')g{^,uy^) (84) 

1 rdv Cdy ^(yy (^^ + l)^^^(^l/w + 1)^^^ 

2 Jo ^ Jo xu,^' u(vv )i/2 

X J'M{l,uy) ^m(^,«1/w) t^:-'^V-''^^'Uuy)g{^,uy^) (85) 
^ dy 1" dy^i'-'yS^:-'y\uy + lY'\uy^ + lf'' 



2 xu^d Jo 

X ^M{^,uy^ J^M^-^.uy^^ -id{uy)G{j-,uy^^ (86) 

/ rfywZ/w(MZ/w+l)'/'j^A/(^,Myw) X (87) 
•J 



X Jo 



where we have defined {y = y/yw) 

2ufidW = r''^ dvv^*-^^'\uvy^ + if/^^M^-^^uvy^) id{uvy^)G{v,uy^) (88) 

To evaluate W we split the integral into f = . . . 1 and v = 1 . . .oo. In the former regime 
we rewrite Q{v^ . . .) in terms of Q{l/v, . . .) using ( l83l) and then transform v 1/v io 
get 

2u^idW = £ dvv^'-'^^/^uvy^ + l)y^T^,(^-^,uvy^) ^diuy„)g{l,uvy^) 

+ J^^^'^ dvv^^-^^/^uvy^ + ifl^TM^^^uvy^) ^d{uvy^)G{v,uy^) (89) 

= j^dvV-^^-^y^ (^ + l)''VM(f^,^)7.(-yw)^(.,^ 

+ f^'""^ dvv^^-^^/^uvy^ + if'^TM^-^.uvy^) id{uvy^)G{v,uy^) (90) 

By performing the w-integrals in fl83|) explicitly, having first inserted the definition fl22l) 
of J-c, one finds for the terms involving Q 

^d{uy^)Q(v,^) _ uy^/v + l 



T^(^d uy^ + 1 



dw w^''-^^/^Tl{w, uyjv)e-^'"^'-'^ (91) 



^(4-d)/2 f'dzdz'izz'Y"-'^/' 

[uy^/v + l){uy^ + 1) Jo 



x(^-^-/ + l)-'^/2(^!^ + l) (^^ + 1) (92) 



Fluctuation-dissipation relations in critical coarsening 



17 



and the replacement u ^ uv gives a similar expression for '^diuvy^)Q{v,uy^) / {T^ad)- 
We now insert these back into fl90|) to obtain 



2^iduW 2(4-'^)/2r(t;/2) 



uy^ + 1 



1 



X dzdz' {zz'f'-^^'^v ~z-z' + l)-^/^ + 1 



uvy 



uy^z' 

V 



+ 1 



X / dzdz'{zz'f''-''^/^{v-z-z' + lY^'\uy^z + l){uy^z' + 1) (93) 

JO 

So far our calculation of C-^^ applies for generic x; to make more progress we consider 
again the large- a; behaviour. In the first f-integral one can use directly the asymptotic 
form ( |55l) of J^m', for the second integral one can show as in the fully magnetized case [I3] 
that the same replacement can be made and the upper integration limit sent to infinity 
thereafter. This gives 



2iidW _d-2 



2 

X (f — 



2(4-d)/2r(^/2)x(2-'^)/4^ 



(d-4)/4 



dvv-''''' / dzdz' {zz'Y^-''^''' 



uy^z 



V 



+ 1 



(94) 



+ 



Then from flSTj) one has 

&'\t,t 



dv / dzdz'(zz'Y''-^^/^(v-z-z' + l)~''/\uy„z + l)(uy„z' + l) 
Jo 

me has 

d-2Y 2(4-^)/2r(ci/2) T>dt 



2 ; mV2(^ + 1)3/2 2fid 



dy^y 



(d-2)/2. 



1 - 



.(d-4)/2 



" POO 


+ J dv . 











(95) 



where the w- integrals are as in flMl) . Carrying out the ?/w-integral, this can be written 



as 



C^'\t,t^ 



d-2 



r((d + 4)/2) u^/H^x-^'^+^y^ 



2 J r((4-rf)/2)r(d + i) (u + 1)3/2 



1 



if we define Vd as in the fully magnetized case, see [Appendix A 

Vd= dv (t;-('^+4)/2 + 1) / dzdz' {zz'f''-''~^'\v -z-z' + l) 



-d/2 



and introduce also the analogous quantities 



V' 



Ad 
rfT2 Ji 



„ 4(rf-l) 



d + 2 



Jo 

dv {v-''/^ + 1) dz dz' izz'Y''~^^/\v -z-z' + l)-''/2 



We have also used in ( l96l) the explicit expression [13] 



-2T, 



2//d (d - 2)(4 - d)2(2-'i)/2r((4 - c/)/2)r((d - 4)/2)r((rf - 2)/2) 



(96) 

(97) 

(98) 
(99) 

(100) 
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With the results (175]) and (jHS]) for the magnetization response and correlation 
in the limit of long, well-separated (x ^ 1) times, we can finally compute the asymptotic 
FDR as 



'd-2 



u + l]{u + iy 



'd{u+ 1) 



r((d + 4)/2) 



X 



UP2(U) 



r((4-d)/2)r(d+l 

'{d+12){u + l) 3u 
2 



4 



^ -m(m + 1)Pi(m) 



(101) 



where we have defined the following 3rd, 2nd and 1st order polynomials in u: 

Psiu) = Au^ + Bu^ + Cu + D 
P2{u) = VdU^ + V:iU + V;i 

Pi(«) = y> + 2\/; (102) 

The general structure of the asymptotic FDR is thus as for (i > 4, i.e. a. ratio of 4th 
order polynomials in u. One can easily check that as — > 4 the coefficients continuously 
approach those for > 4, as they should. Also, the u <C 1-limit of fllOip retrieves the 
prediction for the unmagnetized case in ci < 2 < 4 [5|, [13] 

2 _d-2 

Id ~ d 



(103) 



Conversely, for u ^ 1 one has 



X^ 



d-2 



d + 6 



A + 



' d-2 \ r((d + 4)/2) + 6 

r((4-rf)/2)r(d+i) 4 



v. 



(104) 



which stated in this form agrees with our earlier result for the fully magnetized case [13] . 
(The error was in an incorrect expression for Vd, see Appendix A ) In this regime the 
asymptotic FDR interpolates between X°° = 1/2 for = 2 (as can be shown by using 
that Vd ~ 2/((i — 2) to leading order [T3]) and X°° = 4/5 for d = 4. As is required by 
continuity with the situation for d > 4, the contribution from C^^-* vanishes as c? — 4, 
for any u. An e = 4 — (i-expansion of (110 II) yields 

4(u + l)^ 



X^(u) 



8 + 5u + 30m2 + 20m3 + 

,o 144 + 216m + 48^2 + 160^3 



95u^ + 19u^ 



9(8 + 20m + 30^2 + 20u3 + 5^4)2 ^ 

which in the unmagnetized (m ^ 1) and fully magnetized (m ^ 1) limits reduces to 
X°^{u = 0) = 1/2 - e/8 and X°°{u ^ oo) = 4/5 - (19/450)e, respectively. The 
former value agrees with the well-known result X°° = [d — 2)/d for coarsening in the 
spherical model [5l [13] or the 0{n ^ oo) model [6] from an unmagnetized state. The 
latter, corrected, value now also agrees with the RG calculations for the longitudinal 
fiuctuations of the 0{n oo) model IE 
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One interesting and unexpected feature of fllOip and its expansion fllUSp is that 
the approach to the large-w hmit is non-monotonic: for d close to 4, X°°{u) slightly 
overshoots the limit value "plateau" and then decays down to it, signalling the presence 
of a weak maximum. Expanding (11051) for large u and subtracting off its — oo 
asymptote, one sees that the deviation from the plateau is controlled, to leading order, 
by two terms with opposite signs, scaling respectively as e/u^ and —1/u^. The maximum 
occurs where these two terms compete, that is for u ~ 1/v^, or u = = 0{1). Its 
height above the plateau then scales as e^. To get the scaling function determining the 
shape of the maximum, we therefore normalize the deviation of X°° from the large-w 
plateau by and define 

This scaling function has its maximum at the finite value u = as expected, and 

is positive for u > ^18/5. 

Looking next at dimensions further away from d = 4, Figure [2] shows numerical 
values of for finite u for a few dimensions d between 2 and 4. X°° converges to the 
fully magnetized value (which is near 1/2 for d ^ 2) for large u and to the unmagnetized 
asymptotic FDR {d — 2)/d for u — > as it should. As anticipated, the interpolation 
between these two limits is not, as in > 4, monotonic: initially increases with 
u but "overshoots" its asymptotic limit. This phenomenon becomes more and more 
pronounced as d ^ 2. For d very close to 2, finally, the maximum turns into two poles 
in X°°{u), with X°° being negative in between. 

The first pole is relatively straightforward to analyse from (IIOII) . One needs the 
dependence on 6 = (rf — 2)/2 of ( ITiHTTl) and ( 19714991) for 5 — >• 0. To leading order one 
finds A = 352/4, B = 35/2, C = 1, D = 1/6. The small (5-limits of V^, and V^' one 
gets from (I97H99!) by noticing that the f-integrals become dominated by their large v 
tails as 5 — > 0, giving Vd = 1/6, V/i = 2/6^ and V^' = 1/6^. Gathering these results, 
equation (llOip becomes to leading order 



= {u + lf 



-2u + l)D + h'uV^'^'' .2(^ + 1^3 
' 2 



6--^^ (107) 
2 — u 



which approaches for m — > the unmagnetized limit X = 6 + 0{6^) as it should. For 
larger u we read off that there is a pole at m = 2 beyond which is negative. In fact, 
the expression (I107P shows that X°° ^ —26u'^ for large u whereas we expect convergence 
to the known limit X°° = 1/2. The reason is that the limits u 00 and 5 — do 
not commute: the approach to the eventual asymptotic value takes place on a scale of 
values of u that diverges as 5 — 0. 

The form of the response function as given in (158!) would suggest that the 
appropriate diverging u-scale to consider is m ~ 1/5: in this regime the two terms 
in the square brackets in fl58l) . which cancel exactly for n ^ 00 and d —>■ 2, still give 
a leading order cancellation. However, one finds with a bit of algebra that the limit as 
5 — > of X°^, taken at fixed u' = u6, is simply the constant asymptotic value X°° = 1/2. 
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Figure 2. Asymptotic FDR X°° for tiie magnetization vs u = tyj/Tm, for a 
few dimensions between 2 and 4 as indicated. As d decreases, devefops 
an increasingly pronounced maximum which eventually (sec bottom right graph for 
d = 2.002) turns into two poles separated by a region of negative X°°. 



The crossover to this asymptotic regime must therefore take place on shorter timescales 
u. To explore this, we need to look more closely at the polynomial structure of X°°. As 
observed, X'^ can be written as the ratio of fourth order polynomials, 

^oo ^ g-u^ + bu^ + CM^ + du + e 
a'u^ + b'u^ + c'u'^ + d'u + e' 
The coefRcients can be computed in the limit 5 — > and their leading terms evaluate to 



(109) 



We now consider values of u diverging as some generic power oi 5, u — u"5''^. In the 
limit 5 — > a number of terms can then be dropped: e.g. cu^ in the numerator is always 
subleading compared to bu^ because both b and c are order unity but w ^ 1 for 5 <^ 1. 
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Figure 3. Log-log plot of the absolute value of the asymptotic FDR X°° vs u, for 
6 — 10^"^ (left) and 6 = 10~^ (right). Note that the second pole moves to the right as 
S decreases, with the expected proportionality u ~ S~^^^. Here and in the following 
plots, log = logiQ. 



For the same reason the terms proportional to d, e, c' and e' can never be leading. 
Retaining only the other, potentially leading, terms gives 



U 



/>3 



U" 



(110) 



Comparing powers of 6 shows that the only values of (3 for which in the limit 5 — > 
more than one term survives in either numerator or denominator are (3 = 1/2, 2/3 and 
1. The competing terms at jS = 2/3 are both sub leading so this case is uninteresting. 
Only (3 = 1/2 therefore remains as a non-trivial exponent value to analyse. One then 
has explicitly u" = ub^l'^ and the surviving terms in flllOp are 



25-3/2m"3 _ 1^-3/2 



112 



1 



(111) 



This result matches the magnetized limit = 1/2 for large u!' as it should; for 

= 1/2 it has a pole and for small m" it is negative and small. In the latter regime, 

= — 2ti"^ = —Ibu^ also matches smoothly with the large-w limit of (I107P as it should. 
Figure [3] demonstrates this behaviour by showing on a logarithmic scale the absolute 
value of versus u. As b decreases, the second pole moves to larger m = 1/ (25^/^) as 
expected while the first one occurs at a finite limiting value of m, m = 2. 

It would clearly be desirable to understand in more detail the origins of the highly 
non-trivial behaviour of the asymptotic magnetization FDR for d near 2, and to ascertain 
how this behaviour is reflected in the corresponding FD plots. The response ( l58l) is 
always positive, so from the definition of the FDR in (149|) singularities in X can arise 
only from zeros in dt^C, i.e. from a non-monotonic dependence of the magnetization 
correlator on t„. As we will show, this non-monotonicity arises because the equal-time 
correlator has a pronounced maximum around u = t^/r^ = 1, and this large variance 
of the magnetization fluctuations leaves its imprint in the two-time correlator as a weak 
maximum. 
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The main difficulty we now face is to obtain the behaviour of the correlator also 
for finite x rather than just x ^ 1. This is made possible by the following observation: 
in the S = {d — 2)/2 — > 0-limit, the scaling function J^M{x/y,uy) from (15^ develops a 
non-integrable singularity at y = x. This concentrates the weight of any integrand into 
this region, so that for any function f{y, u) which is smooth dX y = x 

rx ( X \ UX^^"^ 

/ dyJ^M\-^uy\f{y,u)^ —f{x,u) (112) 

JO \y ) ux + 1 

as (5 — 0, i.e. TM{x/y.,uy) acts effectively as J-'M{x/y,uy) = [ux^/'^ /{ux + l)]5(x — y). 
The same observation applies to TuiXlv^i '^Uw), which is obtained by setting x = 1 and 
y = y^. These approximations will yield the leading terms in the correlator in the limit 
5 — > 0. We will also use them for dimensions slightly above 2 to explore numerically 
the x-dependence of the correlator and the resulting FD behaviour. Even though the 
results here no longer have the character of a systematic expansion in 5, they will give 
insights into the non-trivial li-dependence of the FD plots for d close to 2. 

Using the approximation f lll2p . the contributions to C*^^-* coming from 5* vanish 
because the ffist argument of Tuix/y, uy) is always > x > 1. The remaining term F is 
given explicitly in ( 167|) and taking the 5 — limit gives 

C^'Kt,t^) 2x-V^ , (,,^. 

—rj- = 6{ux + imu + 1)3/2 [1 + U5 + U5) (113) 

We note that a naive application of the 5-approximation explained above would in 
this case give an incorrect result, because it produces a leading order cancellation of 
terms when u ~ 1/5. The remaining sub leading term is then of the same order as the 
ffist correction to the ^-approximation. One can nevertheless check that the second 
contribution, 5*, to C^^^ always remains negligible compared to F because it is subject 
to a similar cancellation. 

Next we want to compute C(^). In (IHnD the ffist integral again vanishes because the 
first argument of J^m is always > 1, so 

W=:^^] -, ^ / dz dz' (zz'Y'-'^/' 



f^d \ywj (m?/w + l)(Ma; + 1)3/2 Jo 



X 



— -z-z'+l] {uy^z + 1) {uy^z' + 1) (114) 



which for small 5 evaluates to 

(4-d)/4 ^3-d/2 r 

(My^ + i)(Mx + 1)3/2 r'^' 





In ( ) + i-^^ } (115) 



X J 6'^ X + y 

The i/w-integral in ( l87l) can be performed by again using the 5- approximation and one 
finally gets 

C(2)(t,t^) UX'/^ f 1 , ^ ^ /X + 1 

— ' + 2ud m ' 



Tct 5(u + l)3/2(ux + 1)3/2 la; + 1 ^ ^ 
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+ u'6 



2x2 



a; In 1 -] + In 



X + 1 



(116) 



Adding the results flll3p and (11161) gives the magnetization correlation function for 
generic x and d close to 2. By pulling a factor of x into the curly bracket of the latter, 
one sees that the resulting terms in the bracket are of C(l), 0{6u) and for 
all X, as in (11131) . But the pref actor in (I116p is larger by a factor of m, so for -u ^ 1 
we can neglect C*-^-* against C^'^\ This implies that we can always drop the 0{5u) and 
0{5'^u'^) terms in C*-^-*: either u = 0{1), and then they are subleading compared to the 
0{1) term in C^^\ or m ^ 1 and they are small compared to the corresponding terms 
in C^'^\ Changing then also to a normalization with l/xm instead of 1/t, we can express 
the leading order terms for 5 of the magnetization correlator as a function of the 
scaled times u = t^/r^ and Ut = t/r^^ = xu in the form 



C(t, t„ 
T T 

-"^ c ' m 



11 



U 



5(M+l)3/2(ut + 1)3/2 

,2\ 

+ In 



^In 1-^ 



U + Ut 
U + Ut 



+ 2u6 In 



U + Ut 
Ut 



.Ut - u, 

The equal-time correlator is then obtained by taking the limit u ^ uf. 
C{t,t) ut{4 + Ut + A\n2 6uf + A\n2 6^ul) 



25{ut + lf 



:ii7) 



:ii8) 



Evaluating this numerically for small 5 as shown in Fig. HI we see that it is non-monotonic 
in Ut as anticipated. In fact, the expression (IllSp shows directly that the height of the 
peak at Ut = 0{1) diverges as 1/5 for 5 —> 0, whereas for Ut = 0{l/6) the result 
is of order unity. On the right of Fig. H] we demonstrate that the peak in the equal- 
time correlation function does indeed cause corresponding non-monotonic behaviour 
in the (normalized) two-time correlator C{t,t^) = C{t,ty^)/C{t,t) in the region where 
u = 0{1). Notice that in C the prefactors T^Tj^ that we have isolated on the left of the 
expressions above cancel and we obtain a function of only Ut and u: 

'ut + l\'/' ^u^l-'uT r 2_ 



C{ut,u) 



u + 1 



u 



2x2 



+ u'S 



Ut 



u 



4 + + 41n25u? + 41n252^i3 

,2\ 



U + Ut 



+ 2u5 In 



U + Ut 



Ut 



In 1 



u 



ut 



+ In 



U + Ut 



Ut 



u 



(119) 



We now analyse more closely the nature and scaling of the non-monotonicities of 
the normalized two-time correlator. For Ut of order unity, all terms involving powers 
of 5 can be dropped in (I119p . The resulting function is monotonically increasing in 
u = tw/Tm = . . . M( for > vT3 — 3 ~ 0.61. For larger ut it has a maximum in 
u whose position shifts from vT3 — 3 to an asymptotic limit of 2 as Ut increases. If 
one keeps the terms that are subleading in 5, one sees that C contains a contribution 
scaling as 5'^{ut — u) \n.{ut — u). This yields a term —5^ ln(Mf — u) in the ^-derivative 
of C which diverges to +oo as u ^ Ut, and so gives a positive sign for the derivative 
in the limit. Whenever C has a maximum as a function of u it therefore also has an 
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Figure 4. Left: The equal-time correlator C{t,t) /{TcT^) versus log(Mt) for 6 = 10^^ 
shows pronounced non-monotonic behaviour. Right: Normalized correlator C for the 
same 6 = 10~^, plotted versus log(w) and log(ui). Note the non-monotonicities in the 
M-dependence around log(u) = 0. 




log(u) 



Figure 5. Normalized correlator versus log(M) at fixed ut — 10^ and (5 = 10 ^. 

associated minimum, but this is essentially undetectable as it occurs extremely close to 
Ut, for Ut — u exp(— const/5^). 

Moving to larger Ut of order 1/6, the position of the minimum in C becomes clearly 
separate from Ut- An example of this is shown in Fig. O which graphs the normalized 
correlator as a function of log{u) for a fixed value of Ut with Ut6 = 10: one discerns a small 
maximum followed by a broad minimum. (Numerically, one finds that these features 
merge once d gets sufficiently far above 2, restoring monotonicity.) The maximum in u 
is, for Ut ~ 1/S, always located at m = 2; this matches the behaviour discussed above 
for large Ut of 0{1). Explicitly, if we let 6 tend to in the normalized correlator at fixed 
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Figure 6. Position and scaling of the maximum and the minimum of the normalized 
two-time correlator as a function of u or u" — uS^^^ for a fixed value of ut ~ 1/(5; for 
the example in the plot we have taken Uf6 = ~ 10. 



u and u[ = Ut5 we get 

B(u'u) = limC(u'/6,u) = .,/o.^^^,^^? ; ?k (120) 

The result is shown in Figure [6] (left) and does have a maximum at m = 2 as anticipated. 
This value makes sense since it was also the point where the asymptotic FDR X'^{u) 
diverges for d close to 2. (Given that ~ 1/5 we are automatically in the asymptotic 
regime Ut ^ u.) 

The position Uram of the corresponding minimum of C as a function of u is somewhat 
more subtle. For m[ < 1/2, it is located at Umin ~ 1/5, i-G. Mmin < Ut but with the two 
values being of the same order. As u[ ^ 1/2 from below, Umin^ —>■ 0; for even larger 
values of u[, one finds a different scaling Wmin ~ 5"^^^ so that always Wmin ^ ut. To 
find the minimum position in this regime we need to fix u" = ub^l"^ in the normalized 
correlator when letting 5^0. The typical values of C in this regime turn out to be 
only 0(5^/^) above the plateau B(u[,u — > oo) so we subtract off the latter and divide 
by 5^/^ to define 

M{ut,u ) = hm (121) 

We show a sample plot of this, for a specific value of > 1/2, in Fig. [H] (right). The 
minimum of M occurs at m^j^ = [mJ/(4'uJ — 2)]^/^, which for large v!^ yields m^j^ = 1/2. 
This matches the position of the second divergence of the asymptotic FDR X°°(m), as 
it should. 

As explained above, the derivative of the two-time correlator is always dominated 
by a logarithmically divergent term in the limit u Ut] including prefactors, this 
reads dt„C(t,t^) = —Tc6\n{u[ — u') in the regime Ut = u[/6, u = u'/S. The response 
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function flSSl) is then dominated by the same logarithmic terms: 



R(u'u') = -^ l-u:infl-^) ^ -6\n(u[-u') 
u't [ V ^tJi 



(123) 



The last approximation, which gives the dominant term for u' u[, shows that the 
initial (negative) slope of the FD plot is always exactly equal to one. However, as u[ 
becomes small this becomes undetectable because so does the logarithmic singularity in 
the correlator. 

We are now in a position to analyse the magnetization FD plots for d near 2. 
We consider normalized plots (x vs C) as in d > 4, holding Ut fixed for each plot 
as before to get a valid connection with the FDR X and varying u. We obtain x 
by numerical integration of fl58l) . according to the definition (1501) . and then dividing 
by (11181) . The asymptotic FDR X°°(m) that we have calculated applies in the limit 
Ut ^ u, corresponding to the region in the top left hand corner of an FD plot. Starting 
from the top left corner {u 0) we then expect to see in the FD plots the slopes varying 
as given by X°°{u): initially small (of 0{S)) and negative as usual, then turning positive 
and of order unity, and finally negative again. This S-shape should be present for large 
Ut] for smaller Ut, only part of this variation will be accessible because u < Ut- 

For d close to 2 the above expectations are indeed borne out by numerical evaluation 
as illustrated in Figure [7] for 6 = 10~^. As Ut increases, we start from the fully linear 
FD plot of the unmagnetized case, with negative slope X = X°°{u = 0) ~ 5. A section 
of much larger X then grows and eventually "flips" to the right, producing a region of 
negative FDRs. For much larger values {ut ~ 6^^) the beginning of the FD plot (equal 
times, where the plot meets the horizontal axis) eventually swings back to the left to 
return to the conventional negative slope. The initial slope of —1 also becomes visible. 
At this stage the expected S-shape is complete; it then shrinks progressively towards 
the top left corner as Ut grows and the rest of the plot approaches the close-to-linear 
shape for the fully magnetized case. The region of the plot occupied by the "S" 
scales as l/wf = l/{utSY for ut ^ 1/6. This is clear from (11201) which gives the typical 
values of C at the maximum and minimum, i.e. at the right and left boundary of the 
"S" . The S-shaped region ends where the plot meets the y-axis with an asymptotic slope 
that is Mr independent: this point corresponds to m ^ 0, so we are always in the regime 
Ut ^ u where the asymptotic FDR X°° applies and the negative slope is X°°{u = 0) ~ 5. 
(In fact, this argument applies for any ut, whether or not an actual S-shape is present.) 

The crossover between unmagnetized and fully magnetized behaviour can also be 
seen from the Mr dependence of the |/-axis intercept of the FD plot, which can be 
thought of as its "axis ratio" Y. This is found from the large x limit at fixed Ut of 
the susceptibility, multiplied by and normalized by the equal-time correlator C{t,t) 
from (IllSp . The former is determined from the response function (l58i) by integration, 
x(t,t„) = //^ dt' R(t,t'). Rescaling t' = zt gives 



{ut + 1)1/2 J,^, 




X 



dzz^^-'^/^Utz + lf/^ 1 



Ut + l 



(l_^)(<i-2)/2 



(124) 
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Figure 7. Normalized FD plots for d close to 2, showing normalized susceptibility x 
versus normalized correlation C, for 5 = 10~^ and increasing values of ut as shown 
in each plot. Once the S-shape appears, it remains present for all larger ut but gets 
squashed into a region in the top left hand corner scaling as l/{utS)'^. The slope where 
the plot meets the y-axis is always given by X^{u — 0) S. 
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Figure 8. Axis ratio Y of FD plot, given by the asymptotic normalized susceptibility 
x{ut,u — > 0), versus log(Mt) for S — 10~^. The scaling functions (|125p and (I126P are 
shown and match in the crossover regime (1 ^ <C 1/5) as expected. The dotted 
lines represent their continuations towards larger and smaller ut, respectively. For 
Ut = 0{1), Y is 0{S) as expected from the unmagnetized case, whereas ior ut ^ l/S 
we retrieve the magnetized limit Y — 1/(2 In 2). 



VrTTTt + u;^^^ Iniy/TTt + y/mTt)] (125) 



In the small 5-limit at fixed Ut = 0{1), the square bracket simplifies to l/{ut + 1) and 
the integral can be done explicitly. Multiplying by Tc/C{t,t) gives for the axis ratio 
in this regime 

4 + Mi 

This is of order 6 as expected from the FD plots in Fig. [71 For ut —>■ one gets Y = 6 
exactly, consistent with the known results for the unmagnetized case [13]; for large Ut, 
on the other hand, Y = 26ut. 

In the regime Ut ~ 1/S one finds similarly, by setting ut = u[/5 and taking 5 — 

Y = + (126) 



This is of order unity, again consistent with the FD plots shown above. For -uj ^ 1 it 
approaches 2u[, matching the result from the previous regime, while for u[ ^ oo one 
retrieves F = 1/(2 In 2) in agreement with the result for the fully magnetized case [13j . 
We show the two scaling functions together in Fig. [S], for the example 5 = 10~^. As 
expected the two functions agree in the intermediate regime 1 -C -C 1/5, where the 
axis ratio crosses over from values typical of unmagnetized coarsening {Y ~ 6) to the 
values of order unity for the magnetized scenario. 



5. Discussion 



We have used exact calculations to study the crossover from unmagnetized to magnetized 
initial conditions in the critical coarsening of the spherical ferromagnet. Our focus was 
on the correlation and response functions of the overall magnetization and the associated 
non-equilibrium fiuctuation-dissipation (FD) relations. We derived, in particular, the 
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first exact results (in the nontrivial regime d < 4) for tlie crossover function X°°{u) 
governing the behaviour of the asymptotic FD ratio X°°] u = t^/r^ is the appropriate 
scahng variable, namely the ratio of the earlier measurement time t^ and the timescale 
Tin ~ l/'"^o t)y the initial magnetization mg. While X°°(-u) does interpolate between 
the known unmagnetized {u <C 1) and fully magnetized {u ^ 1) limits, we found that 
unexpectedly the behaviour for intermediate u is not monotonic. In fact, for dimensions 
d ^ 2 close to the lower critical dimension these non-monotonicities turn into pole 
singularities in X°°(m). 

We traced this unusual behaviour to a non-monotonic dependence on the earlier 
time tw of the two-time magnetization correlator C(t, tw), which displays a weak 
maximum at u = 0{1) and a corresponding mininum at u = (for sufficiently 

large Ut = t/Tm = 0{S~^)). We interpreted the maximum as the result of an unusually 
large variance of the magnetization fluctuations in this region, corresponding to a strong 
peak in the equal-time correlator C{t^, t^). The maximum and minimum of C{t, t„) also 
manifest themselves as S-shapes in the magnetization FD plots for d ^ 2. 

As an aside, we notice that non-monotonicities in the asymptotic FDR have 
previously been observed as a function of the lengthscale being probed [l^. Also here 
the effect gets stronger as d 2. On the other hand, the non-monotonic dependence on 
the lengthscale disappears for large enough u at any d > 2, so that it is unclear whether 
the physical mechanism at work here is related to the one causing the complicated 
dependence of on scaled system age that we saw above. 

The calculations for the more general crossover case revealed a typographical error 
in our earlier study of fully magnetized initial conditions [13]. Having corrected this, 
the expansion to flrst order in 4 — d of the asymptotic FDR X°° for the magnetized 
case now agrees with the result of an RG calculation for the 0{n oo) model |14j . 
This agreement suggests, non-trivially, that the spherical and 0{n —>■ oo) models are 
closely related even beyond the leading order Gaussian description of their dynamics. 
One might then suspect similar agreement also with the n-vector model; to verify this, 
it would be desirable to extend existing RG expansions around d = 2 [18] beyond the 
leading term X°° = 1/2 + 0{d — 2). Note that in comparing the spherical with the 
0{n) and n-vector models one has to look at the longitudinal degrees of freedom in the 
latter since these are the ones which — like the magnetization in the spherical case — have 
nonzero average. The transverse fluctuations in coarsening from a state with flnite initial 
magnetization behave differently, giving in the 0{n —>■ oo) model an asymptotic FDR of 
X°° = d/{d + 2) [15j. One expects that the transverse fluctuations in the n-vector model 
would give the same FDR for n — oo. This is consistent with the flrst-order expansion 
X°° = 1/2 + (c? — 2)/8 calculated in [18]. Intriguingly, even though the spherical model 
with its single degree of freedom per lattice site has no direct analogue of transverse 
fluctuations, it gives the same FDR X°° = d/{d -\- 2) for short-range observables when 
the system coarsens from an initially magnetized state. 

In future work, an issue of obvious interest would be to understand how generic 
our results are, i.e. whether similar non-monotoniticies appear also in true short-range 
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models. Field-theoretic calculations near d = 4 |10| for e.g. the 0{n) model should 
in principle be possible, and could be directly compared to the expansion fllOSp of our 
results near d = 4. Our analysis also suggests that if similar expansions were carried 
out near d = 2 [18], very rich behaviour could result. 



Appendix A. Corrections to Ref. [13j 

In this appendix we list the required correction to the relevant equations of [13]. The 
source of the error was equation (8.94) of [13j: it should be replaced by 

r Jdww('^-^y'J^U^)e-'^^-'^- r ^ , 
^^^^ = Git/t^), Gix) = { xJdww('^-^y^:FUw) (A.l) 
^'-y^'^) [ x^''-^y^g{l/x) forx<l 

The old version had an erroneous rather than x^'^~^^^'^ in the second line of the 

curly bracket. All other mistakes are due to trivial propagation of the one above. This 
affects the first integral in each of equations (8.99) and (8.100), whose correct versions 
are 

^= ['duJ^M{x/uy„)u^''-^^/^g{l/u)+ r'""duJ'M{x/uy^)u^''-''^'^G{u) (A.2) 

7d "'0 Jl 

duTM[xuly^)u-^^^'''^I^G{u)^ \ duJ^M{xluy^)u^^-''^'^g{u) (A.3) 



This leads to 
T 

U 



r(^)r(^) 

rx/y. 



[^duJ^M{^) ^'^'^''^''l dy dy'{yy'r-'^/'{l-y-y' + u)-'/' 



+ ^^'-M4^)^''-''^'J/yJ/y'■■■ 



in place of equation (8.104). Equations (8.105,107,108,113) as written in terms of are 
correct, but Vd itself as stated in (8.106) is incorrect. The correct version is ( 1971) in the 
main text. The limiting value V4 for — 4 can be worked out explicitly as V4 = 5/6 
and must replace equation (8.115). This enters X°° (denoted in [13]) only at the 
first order in an e = 4 — ci-expansion 

4 1 9 

X- = --— e + 0(e2) (A.5) 

5 450 ^ ^ ^ ' 

which needs to replace equation (8.116). 

In the opposite limit d ^ 2, the error affects again only subleading contributions, 
so e.g. equation (8.117) stands as written. The first correction oq in the 6 = [d — 2)/2- 
expansion of = 1/6 + Qq + . . . can be obtained as ao = —3/2 (rather than 
Oq = —1/2 — 7r^/12 p^) following the reasoning in [13]. The correct expansion of 
X°° near d = 2 then becomes 

X- = l + ^{d-2) + 0{{d-2r) (A.6) 
which should replace equation (8.118). 



Fluctuation- dissipation relations in critical coarsening 



31 




d 

Figure Al. Asymptotic FDR for the magnetization, for critical coarsening with 
nonzero initial magnetization. Dashed line: old, incorrect version from jl3j : solid 
line: correct version, dotted lines indicate the (corrected) first-order expansions (jA.Sp 
and (|A.6p near d — 4 and d — 2, respectively. 




Figure A2. Normalized magnetization FD plot showing normalized susceptibility x 
versus normalized correlation C in the limit of long times. For d = 2 (lower dotted 
line) and d = 4 (upper dotted line) the old and the corrected versions coincide. The 
deviations between the two versions (correct: full line, old: dashed line) are largest in 
d = 3. The correct plot here is somewhat closer to the straight line obtained for d = 4, 
lying very slightly above it in the right-hand part of the plot. 

Fig. lAll shows the correct d-dependence of compared to the erroneous version 
from [13]. As expected from the discussion above, the quantitative corrections are largest 
around d = 3 and vanish as d approaches 2 or 4. For the sake of comparison, we also 
show in Fig. IA2I the corrected magnetization FD plots: in d = 2 and d = 4 these are as 
before, whereas for d = 3 small quantitative differences are just about visible. 
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Appendix B. Comparison of L^^^ with previous results 

We show briefly in this appendix that the general and exact solution for L^^^ 



dtN{t, t. 



9{U 



9'{t)9{U 



-N,^{t - t. 



(B.l) 



g{t) ' g\t) 

reproduces the long-time results obtained for the unmagnetized and fully magnetized 
limits in [13]. Beginning with c? < 4, because g'{t)/g{t) ~ 1/t the second term on the 
r.h.s. is non-negligible only in the aging regime (t/t^ = x > 1, hence t — tw ~ ^ 1) 



where Neq{t 
L^^\t,t^) = L^S(t 



[2{t-t^)/{A-d)]Lf){t-t 



t^^x 



u + 1 ^ {aux — K,){u + 1) 



from (l40l) . Inserting (fT6!) then gives 
2 x-l 



A-d 



X 



(B.2) 



ux + 1 ' {ux + 1)2 

This result is of the general scaling form L^'^\t,t^) = L''^(t — t^)J^L{x^u)^ with the 
scaling function J^l providing a multiplicative aging correction of the equilibrium 
result. For m ^ 1 or n ^ 1 the M-dependence drops out and we obtain the 
scaling functions found previously for d < A. Speciflcally, for m <^ 1 the square 
bracket simplifles to 1/x and one gets the unmagnetized result J-'iix) = x''^"'^^/^ as 
in [13], whereas for u ^ 1 one retrieves the expression for the fully magnetized case, 
J^l(x) = a;(2-^)/2[2 + (2 - d)x]/{4: - d), also derived in [13]. 

For d > A the situation is a little more complicated because for long times, from fHOj) . 

iVeq(t - tw) = l/f^d + [2(t - tw)/(rf - 4)]L(J(t - t„) (B.3) 

has a constant part of order unity. This means that e.g. in the unmagnetized case, 
where g{t) approaches a constant and one would normally drop the term proportional 
to g'{t) in ( H3l) . a subleading contribution needs to be retained in g'{t). The form of this 
can be found from the Laplace transform iQ together with fj^Tj) : g'{t) has transform of 

^(0) = [s+ 2% - Lg)(o) - (Lg)(.) - mm 



S9[s) 



mo ^ 
s 



m^)/(0) + (l-mg)(/(.)-/(0)) 



Transforming to the time domain gives 

POO 

g'{t)=ml{2T,-Lf^m + ml / dt' Lf^{t') 



;i-^^)/(o)4?w 



(B.4) 



(B.5) 



Here we have neglected the terms arising from /(s) — /(O), which decay as t""'/^ or even 
faster at long times and so will be irrelevant below. We can now systematically analyse 
the order of the various contributions to (IB.ip in the long-time limit, obtained by flxing 
u = t^/Tm and X = t/t^ and taking t^ oo. In the second term of (IB.ip . Nf,q(t — t^) 
scales as 0{tl) + 0{t^^-'^^/^) from flB^D . For g'{t), we note that in > 4 the flrst bracket 



in flB.51) is equal to /i^^. Using also 



c/t^ = cu/t^ and /(O) = c gives for long 



times 



u 



2ux 
d~A 



(B.6) 
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I.e. g'{t) = 0{ut~^) + 0(42-'i)/2). Integrating w.r.t. t yields g{t) = 0{tl) + 0{t^^-'^y^y, 
the leading order term is given explicitly in f|T3|) . Finally, we have L^{t — t^) ~ (t — 
^^^(2-d)/2 ^ ^(2-d)/2(^^_iY2-d)/2 ^ Q(^^(2-d)/2y g^^g ^^^^^ ^^le unmagnetized case 

■u = is special: both terms of (IB.ip then scale as 0{t^~'^^^'^), with Neq{t — t^) = 1/ fXd, 
g{t) = g{t^) = c/ fXd, g'{t) = —cLf^it) to leading order so that 

r - 1 \ ('^-2)/2i 

Lf^it^t^) = - tw) - = Lg)(t - t„) 1 - (B.7) 

Aging effects appear again via the multiplicative correction in the square brackets, which 
agrees with the result derived in [T3] . 

In the magnetized case, the first term of (IB.ip is still of 0{t^^'^^^'^) and given 
by L^^{t — t„){u + l)/{ux + 1). The second term, on the other hand, has a leading 
0{t-^) contribution of {l/fid){u/t„){u + l)/{ux + l)'^ = {l/fidt)ux{u + l)/{ux + l)'^. The 
subleading terms in N^q, g and g' all give corrections to this of relative order 
which compete with the first term of (IB.ip . The overall result can be written in the 
form 

L^'\t,t^) = Lfi{t-t^):Fdx,u) + ^^^^ (B.8) 

In the aging regime (x > 1) the second term dominates; for m — >■ oo it reduces to 
l/{fidtx) = tyj/i^iidt^) consistent with the result of [13]. The full expression for the aging 
correction factor J-'l{x,u) in the first term is rather long so we omit it. At any rate, 
one sees that this first term becomes subleading compared to the second one already for 
time differences t — t„ ~ [t^{u + <^ t^ where J^l{x,u) = = 1. The 

detailed form of the multiplicative aging correction therefore never becomes relevant. 
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